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ABSTRACT 


A  theory  is  developed  for  the  steady  aerodynamic  loading 
on  a  ducted  propeller  with  finite  blade  nuiriber  in  static  and 
low-speed  flight.  The  interaction  of  the  higher  harmonics 
with  the  steady  loading  is  clarified  and  it  is  shown  that,  in 
a  rigorous  sense,  a  high  blade-number  restriction  is  required. 
For  a  shroud  of  small  camber  and  zero  thichness  an  analytic 
soluti.on  is  obtained,  together  with  simple  expressions  for 
shroud  thrust  and  surface  pressures,  stibject  to  a  constant- 
pitch  assumption  on  the  trailing  propeller  vortices.  To  be 
more  precise,  the  pitch  hypothesis  is  dropped  and  an  alterna¬ 
tive  iterative  solution  pursued.  The  final  results,  however, 
of  an  illustrative  calculation  indicate  that  the  two  solutions 
do  not  differ  significantly.  Results  of  a  limited  comparison 
with  experiment  are  quite  reasonable. 
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THE  DUCTED  PKOPELLER 
IN  STATIC  AND  LOW-SPEED  FLIGHT 


Introduction 


In  previous  studies  we  have  treated  the  ducted  propeller 

quite  extensively  for  the  case  of  forward  flight  at  zero 
1-4 

incidence  .  We  now  turn  our  attention  to  the  case  of  static 
and  low-speed  flight,  and  define  the  problem  as  follows:  Given 
the  shroud  geometry  and  operating  conditions,  determine  the 
restilting  thrust  and  time-average  surface  pressures  on  the 
shroud . 

In  general,  study  of  the  static  and  low-speed  problem  has 
been  confined  to  one-dimensional  momentum  analyses,  except  for 
the  work  of  A.  R.  Kri^el^'^  and  H.  R.  Chaplin"^.  Kriebel  has 
obtained  an  expression  for  the  division  of  thrust  between  the 
shroud  and  propeller,  whereas  Chaplin  has  developed  a  nimieri- 
cal  program  directed  mainly  at  computing  the  slipstream  con¬ 
traction.  In  both  analyses  certain  restrictions  are  imposed 
a  priori  which  limit  physical  insight  into  the  complicated 
interactions  that  arise  in  static  and  low-speed  operation. 
Although  we  will  impose  essentially  the  same  restrictions  ex¬ 
cept  for  zero  tip  clearance,  the  nature  of  these  interactions 
will  be  clcirified  in  the  present  development.  More  important, 
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we  will  obtain  convenient  analytic  formulas  for  the  shroud 
thrust  and  s\irface  pressures. 

Analogous  to  o\ir  treatment  of  the  foiward-flight  condi- 
tion,  we  consider  first  the  camber  problem,  using  a  classical 
three-dimensional  vortex  approach.  That  is,  the  shroud  and 
propeller  are  represented  by  a  distribution  of  bound  ring  and 
line  vortices  respectively,  together  with  their  appropriate 
trailing  vortex  systems.  The  effects  of  fluid  viscosity  and 
compressibility  are  neglected  as  well  as  the  presence  of  a 
hi2b,  guide  vanes  and  support  struts.  Once  the  strength  of  the 
bo\ind  shroud  vortices  is  determined,  the  pressures  and  thrust 
may  be  readily  computed. 

It  is  anticipated  that  the  contents  of  this  report  will 
be  condensed  cind  presented  in  a  tabular  format  as  a  simple 
addition  to  the  procedure  of  Ref.  4  . 


CHAPTER  ONE 


FORMULATION 


1.1  Problem  Statement 

We  consider  a  ducted  propeller  in  static  or  steady  low- 
speed  flight  U  at  constant  rotational  propeller  speed  Q 
and  zero  incidence  in  an  unbounded,  inviscid,  incompressible 
fluid  otherwise  at  rest.  The  term  "low-speed"  will  be  made 
precise  later  in  the  sinalysis. 

The  shroud  is  taken  to  be  infinitely  thin  and  of  small 
camber.  The  N  propeller  blades  are  assumed  to  be  infinitely 
slender,  and  the  hiib  radius  is  considered  negligible. 

In  propeller-fixed  coordinates  K,r,0  the  incoming  free 
stream  has  a  uniform  translation  U  in  the  positive  x-direction 
and  a  rotational  component  nr  in  the  positive  0 -direction, 
see  Fig.  1.1  .  The  blades  are  at  rest  in  this  system  but  the 
duct  is  rotating  about  the  propeller  axis  with  angular  speed 
n  .  This  rotation  may,  however,  be  set  equal  to  zero  without 
altering  the  potential  flow  since  the  duct  is  axi symmetric . 

The  resulting  flow  field  is  time -independent  but  does  retain 
a  0 -dependence  because  of  the  finite  blade  number.  Taking 
the  propeller  blade  loading  as  known,  we  seek  to  determine 
the  zeroth  0 -harmonic  of  the  shroud  loading. 
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1.2  Boiandary  Conditions  and  Vortex  Model 

The  flow  field  must  necessarily  satisfy  the  following 
tangent-flow  condition  on  the  shroud, 


tan  e 


v(x,rg,0) 

U+u(x,rg,0) 


(1.1) 


where  €(x)  and  r  (x)  denote  the  local  inclination  angle 

s 

and  radius  of  the  shroud  respectively.  The  components  of 
the  disturbance  velocity  induced  by  the  propeller  and  duct 
u,v,w  are  defined  in  terms  of  the  total  velocity  q  by 


q  =  (u+u)l^  +  vl^  +  (nr4w)ig 


(1.2) 


Where  ij^  #  i^.  ^  i^  are  axial,  radial,  and  circiimferential 
unit  vectors  respectively.  With  R  defined  as  the  duct 
radius  at  the  propeller  plane,  it  follows  from  our  assxmption 
of  small  camber  that 


tan  €  =  0(f:) 

q(x,rg,0)  =  q(x,R,e)(l +0(€)]  (1.3) 

In  addition,  the  magnitude  of  u  lies  between  zero  for  the 
static  case  and  0(u,v,w)  for  low-speed  flight.  The  linear¬ 
ized  form  of  the  tangent-flow  condition  may  therefore  be 
esjpressed  sdLmply  as 
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0  =  v(x,K,e)  (1.4) 

To  compute  the  disturbance  flow,  in  particular  the  radial 
wash  v(x,R,0)  ,  we  introduce  a  classical  vortex  model  to 
represent  the  shroud  and  propeller  as  in  the  forward-flight 
case^,  see  Fig.  1.1  .  That  is,  the  shroud  is  replaced  by  a 
distribution  of  ring  vortices  of  strength  7(x^,6^)  per  unit 
axial  length  over  the  "mean"  duct  r^  =  R  ,  -  c/2  5.  c/2  , 

where  the  v  sxibscripts  will  be  used  in  referring  to  vortex 
elements.  Conservation  of  vorticity  requires  that  these  botind 
vortices  be  accompanied  by  trailing  vortices  of  elemental 
strength  -  ^7/^0^  .  The  propeller  is  represented  by  N 
bo\md  radial  vortex  spihes  of  strength  r(r^)  over  the  blade 
radius,  0  r  £  R  ,  together  with  their  assoc  ated  trailing 

vortices  of  strength  -  r'(r^)  per  unit  radial  length. 
Throughout,  the  prime  will  denote  differentiation  with  respect 
to  the  indicated  argument. 

1.3  Evaluation  of  the  Radial  Wash 

The  radial  wash  v(x,R,0)  contains  0-harmonics  of  order 
0,  N,  2N,  ...  ,  each  of  which  must  satisfy  Eq.  (1.4).  In 
seeking  the  zeroth  harmonic  of 
us  consider  the  zeroth  harmonic  of  Eq.  (1.4), 

0  =  Vq(x,R) 


(1.5) 
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we  may  express  Vq  as  the  siam  of  four  components:  , 

V  '  ,  Vp  and  Vp^  /  which  represent  the  contributions  from 
the  bound  and  trailing  shroud  and  propeller  vortices  re¬ 
spectively.  These  components  can  be  expressed  in  terms  of 
the  unknown  shroud  loading  7  and  the  given  blade  loading 
r  by  means  of  the  Biot-Savart  law. 

This  general  procedure  has  been  carried  out  in  Chapter 

Two  of  Ref.  1  .  The  components  v  and  Vp  may  be  found, 

^0  ^  0 

with  slight  changes  in  nomenclature,  as  the  zeroth  harmonics 


of  and  respectively.  We  fi 


find  that 


(x,R)  = 


c/if 

j  (x-x^)Q^(cJi)  yoi\) 


-c/2 


Vp  (x,R)  =  0 
0 


(1.7) 


where 


=  1  +  (x-x^)^/2R^ 


(1.8) 


and  Qj^  is  the  Legendre  function  of  second  kind  and  one- 
half  order.  The  Legendre  functions  of  second  kind  and  general 
half-integer  order  have  been  tabulated  in  Ref.  8  . 

In  computing  v  ^  we  see  from  the  harmonic  nature  of 
the  shroud  trailing  vortex  system  that 
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V  =  0  (1.9) 

In  order  to  proceed  with  the  evaluation  of  Vp'  ,  it  is 

^  0 

necessary  to  know  the  geometry  of  the  trailing  propeller  vor¬ 
tices.  For  any  given  blade  this  is  determined  by  the  tra¬ 
jectories  of  the  fluid  particles  which  pass  over  the  blade 
and  are  convected  downstream  with  the  total  local  fluid 
velocity  q  ,  given  by  Eq.  (1.2).  In  the  case  of  forward 
flight,  the  perturbational  terms  are  negligible  compared  to  u 
and  Hr  so  that  within  the  confines  of  linear  theory  the 
trailing  vortices  were  assxamed  to  be  regular  helices  with 
pitch  determined  by  U  and  fir^  .  In  low-speed  flight  we  may 
again  neglect  v  and  w  .  For  w  this  is  permissible  because 
it  is  still  quite  small  compared  to  fir  .  For  v  the  justi¬ 
fication  follows  from  the  twofold  action  of  the  shroud:  First, 
by  contributing  an  appreciable  amount  of  thrust  the  shroud 
effectively  reduces  the  required  propeller  disk  loading,  and 
hence  the  resulting  slipstream  contraction,  second  find  more 
important,  slipstream  contraction  is  inhibited  by  the  shroud 
tangent-flow  condition,  particularly  over  the  portion  immedi¬ 
ately  downstream  which,  in  fact,  contributes  most  heavily  to 

Vp'  on  the  shroud.  These  arguments  are  substantiated  by 

°  7 

recent  numerical  results  of  H.  R.  Chaplin'  who  computes  ulti¬ 
mate  contraction  ratios  of  only  a  few  percent. 

Although  we  neglect  v  and  w  ,  we  may  not  neglect  u  , 
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which  is  of  the  same  order  as  u  for  the  low-speed  flight 
condition.  The  axial  component  actually  varies  both  axially 
and  radially,  resulting  in  deformed  helices  of  non-constant 
pitch.  This  causes  considerable  complication.  For  simpli¬ 
city,  we  will  disregard  any  axial  variation  and  approximate 
each  trailing  vortex  as  a  regular  helix  of  pitch  j(r)  .  The 
appropriateness  of  this  important  assumption  will  be  examined 
in  some  detail  later  and,  for  the  present,  we  merely  give  the 
following  motivation:  For  a  free  propeller,  u  increases 
monotonely  from  its  value  at  the  propeller  plane  to  approxi¬ 
mately  twice  that  value  in  the  ultimate  wake,  so  that  a  con¬ 
stant-pitch  assumption  would  clearly  be  inappropriate  in  the 
low-speed  range.  In  the  case  of  the  ducted  propeller,  however, 
the  duct  induces  an  axial  component  which  varies  in  a  comple¬ 
mentary  fashion;  that  is,  it  decays  from  its  value  at  the 
propeller  plane  to  zero  in  the  ultimate  wake.  As  such  it  can 
be  expected  to  remove  much  of  the  axial  variation  in  pitch  so 
that  the  pitch  may,  indeed,  be  sufficiently  constant. 

With  these  approximations,  the  appropriate  expression 

for  on  the  shroud  is  given  by  Eq,  (l6)  of  Ref.  9/  with 

0 

suitable  modification  of  nomenclature,  as 


Vp'(x,R)  = 


N 


0 


47r^R. 


Vr 


A  5^(0^ 

J  J(r  )  '  v'  V 


(1.10) 


where 


Og  =  1  +  [{x-Xp)^  +  (R-r^)^)/2Rr^ 


(1.11) 


and  j  is  the  local  pitch,  as  yet  iindeterroined . 


1.4  The  Resulting  Integral  Equation 

Inserting  Eqs.  (1.6),  (1*7)/  (1‘9)  and  (I.IO)  into  the 
tangent-flow  condition  Eg.  (1.5)  resialts  in  the  following 
integral  equation  in  non-dimensional  form  for  7^  , 


—  r  ^ 

'  47r^  J 


- r  (r  )  dr 


(?v) 


V'  V 


A 


With  fiE  as  oxir  reference  velocity  instead  of  u  ,  which 
P 

may  be  zero,  we  have  non-dimensionalized  according  to 


X  =  x/r 


f  =  r/R 


V  >o/®p 
r  =  r/sjR? 


(1.13) 


and  have  introduced  the  geometric  parameters 


Ul  s  R  /R 

P 

X  =  c/2R 
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(1.14) 


comparing  Eq.  (1.12)  with  the  corresponding  forward- 
flight  equation,  see  Eqs.  (1.4)  and  (1.5)  of  Ref.  2,  we 
observe  a  striding  difference: 
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For  a  given  duct  of  small  geometric  camber,  the 
function  €(x)  appears  explicitly  and  plays  a 
domincint  part  in  the  for^vard-f light  equation, 
but  disappears  completely  from  the  equivalent 
low-speed  equation. 


For  practical  purposes  this  offers  a  tremendous  simpli¬ 
fication  since  the  solution  will  be  independent  of  e (x) 
and  will  depend  only  on  discrete  parameters. 


CHAPTER  TWO 


SOLUTION  OF  THE  INTEGRAL  EQUATION 


2.1  High  Blade-Nurnber  Limitation 


Despite  the  absence  of  the  geometric  camber,  we  point  out 
that  if  j  is  regarded  as  a  known  fimction  then  Eq.  (1.12)  is 
of  the  same  mathematical  form  as  the  forward -flight  equation 

p 

which  has  already  been  inverted  in  closed  form  . 

But  how  is  j  determined?  We  could,  for  example,  take 
j  equal  to  its  actual  value  at  the  propeller  plane,  say  . 
Unfortunately,  the  latter  depends  in  part  upon  and  u^  ^  . 

These  velocities  are  functionals  of  7^  so  that  the  integral 
equation  is  nonlinear.  Since  they  are  f\mctionals  of  the 
higher  harmonics  7^  as  well,  the  integral  equation  is  also 
coupled  with  the  higher  harmonic  integral  equations.  This 
situation  appears  to  be  intractable. 

Instead,  suppose  we  take  j  equal  to  its  actual  value  in 
the  ultimate  wake,  say  .  This  value  still  depends  upon  < 

whereas  u^  tends  to  zero  at  infinity  and  is  not  involved. 

The  velocity  1:1^,  is  a  complicated  functional  of  the  higher 
harmonics  7^  ,  but  is  independent  of  7^  so  that  the  result¬ 
ing  integral  equation  is  linear.  However,  the  coupling  re¬ 
mains.  It  appears  that  the  only  hope  of  arriving  at  a  tract¬ 
able  situation  is  in  further  assuming  the  blade  number  to  be 
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infinite,  thus  forcing  ,  and  hence  the  coupling,  to  zero. 

This  restriction  is  imposed  for  the  sa3ce  of  mathematical 
rigor;  practically,  we  would  e^qoect  the  results  to  be  quite 
valid  for  blade  numbers  as  low  as  3  4  . 


2.2  Determination  of  Pitch 

To  evaluate  j„(?)  we  note  that  Up  is  identically 
zero,  and  write 


U+Up' 

Hr  =  J  + 
P 


(2.1) 


where  Up^  is  evaluated  at  x  =  «  . 

To  deal  with  the  Up^  term,  we  must  specify  the  blade  cir- 
culation  r(r)  .  As  in  the  case  of  forward  flight  ,  we  choose 
a  distribution  such  that  the  Betz  displacement  condition  is 
satisfied  in  the  ultimate  waXe.  For  our  case  of  high  blade  num¬ 
ber  the  Betz  condition  requires  a  uniform  jet  velocity  at  x  =  »  . 
It  follows  that  Up'  in  Eq.  (2.1),  and  hence  ,  must  be  con¬ 
stant.  From  Eq.  (22)  of  Ref.  9  we  have  at  x  =  »  , 


r  >  M- 

f  <  |i 


(2.2) 


and  since  both  Up'  and  are  to  be  independent  of  r  in 

the  ultimate  wake,  we  must  have  r(f)  equal  to  a  constant. 
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say  r  .  This  means  that  all  the  trailing  vorticity  is  concen¬ 
trated  on  the  slipstream  surface  r  =  |i  .  The  pertinent  value 
of  Up'  to  be  used  in  Eq.  (2.1)  is  therefore  the  value  ^  the 
slipstream.  Unfortunately,  this  value  is  not  available  from 
Eq.  (2.2).  Instead,  we  proceed  by  computing  the  propeller 
thrust  T  in  two  different  ways: 

First  define  the  propeller  thrust  coefficient, 

S  T/%p{nEp)^(irR2)  (2.3) 

From  the  Kutta-Joukowski  relation,  consistent  with  our  assump¬ 
tion  that  w  «  flr  ,  we  have 

R 

C.J,  =  p  J  ^  nrNP  dr  /h-npa^R^  =  ^  (2.4) 

0 

On  the  other  hand,  momentiam  considerations  lead  to 

Cj  =  (j  +  5j,-)^- (2.5) 


where  Up'  here  denotes  the  value  inside  the  jet  at  x  =  oo  , 
or  C^2j^  according  to  Eqs.  (2.2)  and  (2.4).  The  resulting 
Eq.  (2.5)  may  be  solved  explicitly  for  and  we  find  that 


[j  +  Jcj  + 


(2.6) 
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Using  this  resxilt  in  Eq.  (2-1),  we  find  as  we  might  suspect 
that  the  appropriate  value  of  Up^  on  the  slipstream  must  be 
the  average  of  the  "inner"  and  "outer"  values,  . 

It  is  convenient,  at  this  point,  to  give  precise  defini¬ 
tions  of  the  so-called  "high-"  and  "low-"  speed  flight  condi¬ 
tions.  Roughly  speaking,  high-speed  flight  (also  called  "for¬ 
ward  flight"  or  "cruise")  corresponds  to  large  j  and  small 
Cy  .  In  the  limiting  case  we  have  from  Eq.  (2.6)  that  -  J  , 
and  so  a  measure  of  our  "closeness"  to  ideal  forward  flight  is 
the  closeness  of  to  unity.  Observing  from  Eq.  (2.6) 

that 

jyj  »  1  +  cy4j2  (2.7) 

we  have  cy4j^  <<  1  as  our  criterion.  Analogously, 

^00  "*  ^  low-speed  limit,  and 

jy(  n^2)  «  1  +  J/V^  (2.8) 

Our  criterion  for  low-speed  flight  may  therefore  be  stated  as 
J/  «  1  .  Both  conditions  may  be  given  in  terms  of  a 
single  "flight  parameter"  6  =  j/ ,  as 

HIGH-SPEED:  1/46^  «  1 


LOW-SPEED: 


6  «  1 


(2.9) 
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2.3  solution  of  the  integral  Equation 

With  f(r)  constant  and  j(r)  given  by  of  Eq.  (2.6), 

the  integral  equation  reduces  to 


s  1  r 


with 


=  1  +  [(x-Xp)^+  (l-n)^]/2p. 


(2.11) 


In  addition,  we  have  dropped  the  zero  subscript  on  7  since 
we  have  forced  the  higher  harmonics  to  zero. 

Following  Ref.  2,  we  introduce  a  new  axial  coordinate  0 
such  that 


X  =  -  X  cos  0 


(2.12) 


and  seek  a  Glauert  series  solution  of  the  form 


ggCOt|  +  /  sin  n't 


(a.i3) 


If  we  express 


00 

Q^(®3)  =  1o  -  I  “Jn 


(2.14) 
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for  convenience,  the  solution  is  given  in  matrix  form  as 


{g}  =  ^  [I-P]“^(q) 

-*co 


(2.15) 


where  (g)  and  (q)  denote  the  infinite  column  vectors 


(g)  =  i  g. 


(q)  n 


(2.16) 


TlJ  being  the  identity  matrix,  with  imity  for  each  element 
on  the  main  diagonal  and  zero  for  every  other  element,  and 
[p]  the  "curvatiire  matrix"  of  Ref.  2,  the  elements  of  which 
depend  upon  X  and  tend  to  zero  in  the  two-dimensional  limit 
of  X-»  0  . 

2 

It  has  been  shown  that  the  Neimann  e^ansion 


=  (ij  +  [p]  +  [p]^  + 


(2.17) 


is  valid  for  all  finite  values  of  X  .  The  inverse  matrix 
has  been  evaluated^  accordingly  for  X  =  0.25/  0.50,  0.75 
and  1.0,  and  is  reproduced  here,  except  for  X  =  1.0,  in 
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Tables  1.1  -  1.3  • 

The  q^(  X,  p.,Xp) ' s  have  been  evaluated  numerically  to 

within  iO.OOl  and  are  tabulated  for  n=0/  1,  2,  ...  12  ; 

x  =  0.25,  0.375,  0.50,  0.625,  0.75  ;  =  0.90,  0.9^,  0.97, 

1.0  ;  and  x  /c  =  0.0,  ±0.125,  ±0-25,  ±0.375,  in  Tables  2.1  - 
P 

2.5  •  The  "in-between"  values,  X  =  0.375  and  0.625,  were 
included  to  retain  a  consistent  degree  of  accuracy  in  inter¬ 
polation,  and  the  value  p.  =  1.0  to  facilitate  interpolation 
for  cases  of  small  tip  clearance. 

In  siammary  then,  the  static  and  low-speed  solution  is 

given  by  Eq.  (2.13)  together  with  Eq.  (2.15)  and  Tables  1.1  - 

2 

2.5  ,  cf.  the  forward-flight  solution  . 


CHAPTER  THREE 


SHROUD  SURFACE  PRESSURES  AND  THRUST 


3.1  Shroud  Surface  Pressures 

We  define  the  shroud  surface  pressure  coefficient,  based 
upon  the  propeller  tip  speed,  as 


_  P-Pq. 


(3.1) 


where  p  and  p^  are  the  static  pressures  on  the  shroud  and 
at  infJLnity  respectively.  From  Bernoulli’s  equation  we  have 

P-P«  =  -  J5P[(U+u)^-U^]  (3.2) 

since  v  »  0  ,  and  Eq.  (3*1)  becomes 

Op  =  -(2J  +  u)3  (3.3) 

We  call  attention  to  the  nonlinear  term  u  .  For  finite 
blade  nvunber  the  higher  harmonics  of  u  would  therefore  con¬ 
tribute  to  the  zeroth  pressure  harmonic  since  terms  of  the 
form  cos  m0  =  (ij  +  ijcos  2me)  and  sin  m0  =  (%-^scos  2m0) 
contain  axisymmetric  portions.  This  coupling  with  the  higher 
harmonics  in  the  zeroth  pressure  harmonic  would  be  in  addition 
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to  the  coupling  of  the  governing  integral  equations,  discussed 
above  in  Section  2.1,  and  once  again  emphasizes  the  need  to 
limit  our  attention  to  high  blade  numbers. 

Now,  on  the  shroud  u  is  equal  to  the  sum  of  and 

Up'  .  By  means  of  a  careful  limiting  process  it  may  be  shown, 
see  Eq.  (I.IO)  of  Eef.  3/  that 

X 

(3.4) 

The  first  term  on  the  right  hand  side  represents  the  well-known 
two-dimensional  contribution,  and  the  second  term  supplies  the 
effect  of  three-dimensionality.  From  Hough  again,  we  may 
express  u^  conveniently  as 

00 

=  +  -^  cot^  ^  X  J  {0)'‘!s](g) 

(3-5) 

where  the  [s(X)]  matrix  is  tabulated  in  Ref.  3  and  repro¬ 
duced  here  in  Tables  3*1  --3.3  fon  X  =  0.25/  0.50,  and  0.75  • 
The  transpose  {0}  vector  is  defined  as 

(0}*^  =  H  (1/  -  COS0,  -  cos20,  —  )  (3.6) 

Turning  to  Up,  ,  we  find,  from  Eq.  (28)  of  Ref.  9, 


2i 


(3.7) 

where  is  Hevcman's  lambda  function  with  argument  3 

and  modulus  k  defined  as 


(3.8) 

This  fxmction  was  tabulated  in  Ref.  10,  and  a  shorter 

tabulation  appears  in  Ref.  11  .  For  p.  and  x  in  the  range 

P 

of  practical  interest  k  is  slightly  less  than  or  equal  to 
unity  over  the  shroud,  and  the  simple  approximation 

A^(3,k)  «  A^(3,l)  =  23/7r  (3.9) 

is  sufficiently  accurate. 

With  u  computed  as  the  sura  of  and  Up'  ,  from 

Eqs.  (3*5)  (3-7)/  the  pressure  coefficient  on  the  shroud 

is  then  given  by  Eq.  (3.3).  In  addition  we  may  compute  the 

net  surface  pressure  coefficient,  c  ]  =  c  (x,!*”)  -  c  (x,!*^)  , 

P  P  P 

as 


(3.10) 

The  first  term  -2J7  coincides  with  the  corresponding 


Up  4^/1)  = 


SttJ 


x-x 
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forward-flight  result,  and  the  additional  two  terras  represent 
the  nonlinear  contributions.  Physically,  Eq.  (3.IO)  expresses 
the  Kutta-Jo\ikowski  interaction  between  the  axial  velocity 
and  the  bound  shroud  circulation. 


3.2  Shroud  Thrust 

From  Lagally's  extension  of  the  general  Kutta-Joukowski 

formula,  the  elemental  thrust  dt  ,  positive  upstream,  on  a 

12 

shroud  ring  vortex  of  strength  is 


dt  =  -  p(vpO(27rR  7dx^) 


(3.11) 


Defining  the  shroud  thrust  coefficient  as 


Ct  5  t/!sp(nRp)®(7rR2) 


(3.12) 


we  find,  with  Vp'  from  the  left  side  of  Eq.  (2.10),  and  7 
and  Qj^  from  Eqs.  (2.13)  (2.l4)  respectively,  that 


Q  p 

^  J  y 

-'00^  -X 


xc 


T 


4j 


00 

-31190  -Z'In(9n+l  '  9n-l) 


n=l 


(3.13) 


q  ‘ s  in  terms  of  the  g  ' s  by 


If  we  express  the 
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means  of  Eq.  (2.15)  and  recall  the  properties  of  the  [p] 
matrix  from  Section  2.2  of  Ref.  2,  we  find  that  reduces 

to  the  single  term 

^  (3.14) 

This  we  recognize  immediately  as  the  classical  leading-edge 
suction  force.  It  is  not  unexpected,  of  course,  since  -y  is 
independent  of  the  geometric  shroud  camber  in  the  present 
theory.  In  other  words,  the  shroud  thrust  is  the  same  as  if 
the  camber  were  identically  zero,  in  which  case  the  only 
possible  axial  force  is  the  leading -edge  suction  force.  This 
result  was  anticipated  several  years  ago  by  A.  H.  Sacks,  in  a 
private  communication,  based  upon  physical  reasoning. 


CHAPTER  FOUR 


INVESTIGATION  OF  CONSTANT-PITCH  HYPOTHESIS 


4.1  Pitch  Iteration 

To  investigate  the  degree  of  validity  of  the  constant- 
pitch  assmnption,  we  will  set  up  an  iterative  solution  which^ 
in  the  limit,  is  expected  to  yield  the  exact  pitch  distribu¬ 
tion  including  the  axial  variation.  At  the  same  time,  it  should 
also  provide  considerable  physical  insight  into  the  problem. 

If  we  take  as  our  initial  value  j  ,  that  is 

C^4j^  and  0  ,  we  find  the  resulting  shroud 

loading  explicitly  from  Eq.  (2.15).  With  this  known, 

we  may  compute  a  new  pitch  distribution  j^^^(x)  according  to 

j^^^(x)  =  J  +  U^^^(X,M-)  +  Ujl^^(x,M^)  ,  £  X  <  CO  (4.1) 

where  u^'‘'(x,|x)  is  taken  as  the  average  of  the  "inner"  and 
"outer"  values  as  before.  Integration  of  over  the  shroud 

yields,  cf.  Eq.  (52)  of  Ref.  9  , 
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The  propeller-induced  contribution  may  be  found  from  Eq.  (28) 
of  Ref.  9  •  After  averaging  and  changing  nomenclatiare,  we  have 


u^'^(x,M.)  = 


<3^5(05) 


(4.4) 


where 


=  1  + 


If  the  right  hand  side  of  Eq.  (4.2)  is  integrated  numeri¬ 
cally,  then  j^^^(x)  is  determined  at  any  x  location.  In 
the  limit  as  x  -*■  «  ,  “*■  ^  and  C,p/4j^  so  that 

-  J„  • 

With  the  new  trailing  vortex  geometry  thus  specified,  the 
next  step  is  to  compute  the  corresponding  propeller-induced 

wash  V^^^  3^”  +•>»<»  aVi-rr>ii/^  Viom-io  +*>>0  noti/  «VirT»iir1  trrsT— 


on  the  shroud,  and  hence  the  new  shroud  vor¬ 


tex  loading  7'  .  Carrying  this  out,  we  find  the  result  is 

equivalent  to  the  integration  over  a  semi-infinite  cylinder  of 
ring  vortices  of  strength  per  unit  length  Nr/27rRpj (x^)  , 


(4-6) 


(x  ) 


For  practical  purposes  we  split  the  integration  into  two  ranges: 

the  first  from  x  to  X  ,  and  the  second  from  X  to  «  .  we 

P  /  N 

choose  X  large  enough  so  that  j'  ' [x)  may  be  approximated 


5 
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by  in  the  latter  "tail"  portion.  The  integration  in 

this  portion  may  therefore  be  performed  analytically,  and 


(•  p  x-x 

P 


J5 

C 


1  + 


(x-X)^+  (l-n)^ 
2H 


('t-7) 


cf.  the  left  hand  side  of  Eq.  (2.10)  with  x^  replaced  by  X  . 

We  may  now  determine  7^^^  .  Recalling  again  the  left 
hand  side  of  Eq.  (2.10),  or 


=  -  ^015(03) 


(4.8) 


as  well  as  Eqs.  (2.14)  cind  (2.15),  we  see  that  if  we  expand 


.  ^  =  q(l)  _  y  q(^)co3  n0 


(4.9) 


n=l 


then  7^^^  is  given  by 


(4.10) 


and  the  Glauer>\  series  of  Eq.  (2.13)  with  g 


(1) 

0 


instead  of 
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and  so  forth. 

This  process  may  be  repeated.  The  only  change  necessary 
is  in  the  evaluation  of  •  Since  the  pitch  is  no  longer 

constant,  Eq.  (4.4)  does  not  apply.  For  and  all  sub¬ 

sequent  iterates,  we  have  the  analogous  form  to  Eq.  (4.6) 
for  the  radial  component,  namely 


u^?^(x,n) 


where 


[04(05) -QIjjCos)!  '^v 


(4.11) 


Og  =  1  + 


(4.12) 


4.2  Illustrative  Results 

We  have  carried  out  tlie  above  process  on  the  CDC  1604  at 
the  Cornell  University  Computing  center  for  the  typical  case 

J  =  0 

Cr^  =  0.1 

X  =  0.5 
\L  =  0.9 

Xp  =  0  (4.13) 

which  will  be  considered  further  in  the  next  chapter.  Here  we 
are  only  concerned  with  the  iteration.  The  static  case  was 
chosen  since  it  is  the  most  critical;  that  is,  the  pitch  hypothe- 


28 


sis  can  be  expected  to  be  least  acciirate  at  the  static  condition 
and  to  become  more  acciirate  as  the  flight  speed  is  increased. 

The  results  are  shown  for  a  two-step  iteration  in  Fig.  4.1  . 
We  see  immediately  that  the  axial  variation  of  and  Up^ 

are  quite  "complementary"^  resulting  in  a  reasonably  constant 
pitch  distribution.  We  also  see  from  a  comparison  of 
and  that,  qualitatively,  the  effect  of  the  axial  varia¬ 

tion  in  pitch  is  to  shift  the  "effective"  propeller  plane  rear¬ 
ward  by  approximately  2^  of  the  chordlength. 

In  judging  the  validity  of  the  constant -pitch  assumption, 
however,  we  are  not  concerned  with  the  discrepancy  between  j(x) 
and  per  se  but  rather  with  the  resulting  discrepancy  in 

the  final  quantities  of  interest,  such  as  the  shroud  loading  7 


and  the  shroud  thrust  C,  . 

t 

n  = 12  ,  the  coefficients  of 


Truncating  the  Glauert  series  at 
7^^^  and  7^^^  are  given  by 


r 

1 

r 

0.1349 

0.1365 

0.0215 

0.0348 

0.0631 

0.0623 

-0.0001 

-0.0065 

-0.0214 

-0.0191 

0.0000 

^  (g(l))  =  < 

0.0019 

<  0.0096 

1  0.0078  !>■ 

0.0000 

-0.0005 

-0.0049 

-0.0039 

1 

0.0000 

0.0002 

0.0026 

1  0.0022 

o.ooc^ 

-0.0001 

-0.0014 

-0.0012 

L 

(4.14) 
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and  and  are  compared  in  Fig.  4.1  .  While  the 

relative  differences  in  some  of  the  coefficients  are  appreciable, 
agrees  with  quite  closely.  Furthermore,  g^^^ 

and  g^  ^  differ  by  only  1^  so  that,  according  to  Eg.  (3.13), 
and  differ  by  only  2^  . 


CHAPTER  FIVE 


NUMERICAL  EXAMPLES 

5.1  Input  Data 

Let  us  consider  the  static  case  J  =  0  with  the  following 
laical  values  for  the  parameters, 

=  0.1 
X  =  0.5 
p  =  0.9 

Xp/c  =  0.0  ,  -0.25  (5-1) 

where  the  first  case,  for  x  /c  =  0.0  ,  is  identical  to  the 

P 

one  considered  above  in  section  4.2  .  The  camber  is  asstoned 
to  be  small,  but  arbitrary,  since  it  has  no  effect  in  the 
present  theory. 

5.2  Results 

In  order  to  determine  the  shroud  surface  pres store  coeffi¬ 
cient  and  the  shroud  thrust,  we  first  compute  the  vortex  loading 
7  .  We  find  (g)  from  Eq.  (2.15)  and  Tables  1.2  and  2.3  • 

For  x  /c  =  0.0  the  result  is  given  by  the  first  of  Eqs.  (4.l4) 
and  is  reproduced  below  for  convenience,  together  with  the 

result  for  x  /c  =  -0.25  / 

P 
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r 

0.1349 

0-1327 

0.0215 

-0.0499 

0.0631 

0.0312 

-0.0001 

j 

0.0317 

-0.0214 

1 

1 

i  1 

0.0087 

0.0000 

1 

1 

1 

!  -O-OO69 

<  0.0096 

>  .  < 

-0.0081  1 

0.0000 

-0.0023 

-0.0049 

0.0023 

0.0000 

1 

0.0027 

0.0026 

0.0008 

0.0000 

-0.0009 

-0.0014  i 

-0.0011 

(5-2) 


respectively.  The  corresponding  >'3  are  computed  from 
Eq.  (2.13)  and  plotted  in  Fig.  5-1  • 

The  inner  and  outer  shroud  surface  pressures  are  then 
calculated  from  Eqs.  (3-3).  (3-5)/  (3-6),  (3-7),  (3-8),  (3-9) 
and  Table  3-2  ,  and  the  results  are  shown  in  Fig.  5-2  .  Two 
striking  features  of  these  curves  are  apparent:  first,  the 
outer  surface  pressure  distribution  is  virtually  independent  of 
the  propeller  position  and,  second,  the  inner  surface  pressure 
coefficient  decays  to  zero  almost  immediately  downstream  of 
the  propeller  plane.  This  lack  of  "communication",  expressed 
by  the  first  result,  is  not  \inexpected,  on  physical  grounds. 
Regarding  the  second  result,  we  find  that  on  the  inner  surface 
Up^  and  iVy  are  additive  upstream  of  the  propeller  plane, 
whereas  they  oppose  each  other  downstream  of  it. 
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consequently,  it  is  natural  to  wonder  whether,  for  practi¬ 
cal  purposes,  we  may  replace  the  two  parameters  X  and  x 

P 

hy  the  single  parameter  (x  +X)  ;  that  is,  the  nondimens ional 

Jb' 

distance  from  the  leading  edge  to  the  propeller  plane.  The 
results  of  sample  calculations,  however,  are  found  to  differ 
significantly  so  that  the  above  simplification  is  not  valid. 

Finally,  from  Eq.  (3.13) »  the  shroud-to-propeller  thrust 
ratio  obtained  for  the  two  cases  of  Eqs.  (4.13)  and  (5-1)  are 


Cf/CT 


0.683 

0.706 


Xp/c  =  -O.25 

X  /c  =  0.0 

P 


(5.3) 


The  closeness  of  these  two  results  seems  to  support  Kriebel's 
a  priori  assumption^ that  C^/C,p  is  independent  of  x^/c  . 
But  we  emphasize  that  this  is  not  the  case,  even  in  an  approxi¬ 
mate  sense,  for  the  detailed  shroud  surface-pressure  distribu¬ 
tion,  see  Fig.  5*2  . 


5.3  Comparison  with  Experiment 

Although  the  present  theory  is  incomplete  in  the  sense 
that  the  effects  of  shroud  thickness  are  not  yet  incorporated, 
it  is  nevertheless  of  interest  to  compare  it  with  experiment. 

In  the  absence  of  suitable  detailed  surface  pressure  measure¬ 
ments,  we  are  limited  to  only  a  gross  comparison  based  upon  the 
shroud  thrust.  If  we  use  the  inner  shroud  radius  at  the  pro- 
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peller  plane,  comparison  with  the  results  of  K.  J.  Grunwald 

and  K.  W.  Goodson^^  where  J  =  0  ,  X  =  0.620  ,  x  /c  =  -0.l46 

P 

and  |i  =  0.995  !  gives  a  predicted  value  of  which  is 

about  30^  lower  than  their  measxired  value,  since  is  a 

second-order  quantity,  though,  the  corresponding  error  in  the 
shroud  circulation  is  only  half  as  great;  i.  e.  about  15^  , 
which  is  quite  reasonable. 


CONCLUSIONS 


A  theory  has  been  developed  for  the  steady  aerodynamic 
loading  on  a  ducted  propeller  with  finite  blade  number  in 
static  and  low-speed  flight.  For  a  shroud  of  small  camber 
and  zero  thickness,  we  have  shown  that: 


The  effect  of  the  geometric  camber  is  of  higher  order 
and  may  be  neglected. 

The  steady  shroud  circulation  is  coupled  nonl inear ly 
with  the  higher  harmonics  through  the  pitch  of  the  pro¬ 
peller  Wcike.  This  coupling  imposes,  in  a  rigorous 
sense,  a  high  blade-number  restriction. 

An  analytic  solution  together  with  simple  expressions 
for  the  surface  pressures  and  shroud  thrust  can  be 
obtained,  subject  to  a  reasonable  assumption  of  constant 
pitch  for  the  trailing  propeller  vortices. 

The  pitch  hypothesis  may  be  dropped  and  an  iterative 
solution  pvirsued,  but  an  illustrative  calculation  indicates 
that  the  corresponding  final  results  do  not  differ 
significantly . 

The  outer  shroud  surface  pressure  is  q\iite  insensitive 
to  the  axial  propeller  position,  and  the  inner  shroud 
surface  pressure  decays  almost  to  zero  immediately 
behind  the  propeller  plane. 

Results  of  a  limited  comparison  with  experiment  are 
quite  reasonable. 


It  is  anticipated  that  the  contents  of  this  report  will 
be  condensed  cuid  presented  in  a  tabular  format  as  a  simple 
addition  to  the  procedure  of  Ref.  4  . 
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table  1.2 

(TlJ  +  [P]  +  [p]2  +  ...  )  MATRIX  FOR  X  =  0-50 


f 


)  KATRJOC  FOR  \  *0.75 


BLANK  PAGE 


Xp/c  =  0.0 


X  /c  =  ±0.125 


n  \. 

0.90 

0.94 

0.97 

\ 

1.00 

0.90 

0.94 

0.97 

\ 

1.00 

0 

1.744 

1*913 

2.044 

2.178 

1.735 

1.907 

2.042 

2.180 

1 

1 

0 

! 

0 

0 

0 

±0.293 

1 

±0.363 

±0.422 

±0.484 

' 

2 

0.440 

1 
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0.769 

0.982 
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0.667 

0.858 

3 

0 

0 

0 

0 
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- 

4 

-0.104 

-0.192 
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-0.104 
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i 

6 
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0.080 
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0.011 
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- 

7 

0 

0 

0 

0 

+0.017 

+0.051 
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+0.280 

± 

8 

-0.011 

-0.037 

-0.096 

-0.250 

0.003 

0.013 

0.039 

0.109 

9 

0 

0 

0 

0 

±0.005 

±0.020 

±0.057 

±0.169 

+ 

10 

0.004 

0.019 

0.060 

0.200 

-0.003 

-0.013 

-0.047 

-0.163 

11 

0 

0 

0 

0 

+0.001 

+0.005 

+0.017 

+0.064 

± 

12 

-0.002 

-0.010 
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0.001 
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0.037 
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1 

— 1 

TABLE  2.1 
COEFFICIENTS  FOP 


X  /c  =  ±0.25 


/  s/ 


0.90 

0.94 

0.97 

1.00 

0.90 

1.705 

1.887 

2.053 

2.185 

1.639 

±0.557 

±0.703 

±0.831 

±0.968 

±0.746 
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0.296 
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-0.006 
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0.125 
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1 
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0.000 
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X  /c  =  ±0.375 


0.94 

0.97 

1.00 
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±1.190 

±1.454 

-0.017 

-0.076 
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0.242 
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0.166 
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0.90 

0.94 

0.97 

\ 

1 

1.00 

✓ 

0.90 

0.94 

0.97 

1.00 

1 

1.462  1.595  1.698  i.8o4 
±0.645  ±0.756  ±0.846  iO.939 
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x^/c  =  ±0.25  Xp/c  =  ±0.375 


0.90 

0.94 

0.97 

1.00 

0.90 

0.94 

0.97 

1.00 

1.277 

1.382 

1.463 

1.545 

1.244 

1.366 

1 .465 

1.568 

±0.677 

±0.766 

±0.835 

±0.906 

±0.944 

±1.100 

±1.230 

±1.367 

0.287 

0.341 

0.393 

0.456 

-0.056 

-0.097 

-0.130 

-0.157 

±0.328 

±0.432 

±0.533 

±0.658 

±0.l8i 

±0.229 

±0.284 

±0.367 

0.091 

0.139 
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0.250 

0.158 

0.241 

0.338 

0.483 

±0.071 

±0.104 
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±0.199 
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±0.224 

±0.356 
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-0.219 

-0.333 
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±0.091 

±0.004 

±0.007 

±0.010 

±0.017 

-0.011 

-0.032 

-0.073 

-0.167 

0.005 

0.017 

0.043 

0.122 

LE  2.3 

ITS  FOR  X  =  0.50 


X  /c  =0.0  X  /c  =  ±0 . 125 

P  P 


X 

0.90 

0.94 

0.97 

\ 

1.00 

y 
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±0.347 

±0.381 
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±0.433 
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0.719 
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0 

0 

0 

0 
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+0.330 

+0.386 

+0.449 
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0 

0 

0 

±0.164 

±0.230 

±0.295 

±0.379 

6 

0.127 

0.186 

0.248 

0.332 

0.009 

0.011 

0.014 

0.018 

7 

0 

0 

0 

0 

+0.088 

+0.140 

+0.197 

+0.279 

8 

-0.069 

-0.116 

-0.170 

-0.249 

0.027 

0.048 

0.073 

0.109 

9 

0 

0 

0 

0 

±0.039 

±0.070 

±0.108 

±0.169 

10 

0.040 

0.076 

0.123 

0.200 

-0.031 

-0.060 

-0.099 

-0.163 

11 

0 

0 

0 

0 

+0.012 

+0.022 

+0.038 

+0.064 

12 

-0.025 

-0.053 

-0.094 

-0.166 

0.023 

0.050 

0.091 

0.166 
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X  /c  =  ±0.25  X  /c  =  ±0.375 

Jr  Jr 


0.90 

0.94 

0.97 

1.00 

0.90 

0.94 

0.97 

l.CO 

1.132 

1.219 

1.286 

1-354 

1.114 

1.218 

1.300 

1.384 

±0.683 

±0.756 

±0.813 

±0-871 

±0 . 966 

±1.099 

±1.207 

±1.319 

0.295 

0.344 

0.388 

0.440 

-0.086 

-0.120 

~0.l46 

-0.167 

+0.371 

+0.465 

+0.551 

+0.653 

+0.197 

1 

+0.244 

+0.293 

+0.363 

0-112 

0.157 

0.199 

0.250 

0.193 

0.274 

C.36I 

0.482 

±0-085 

±0.117 

±0,151 

±0.198 

+0.102 

+0.169 

+0.246 

+0.356 

-0.110 

-0.171 

-0.238 

-0.332 

0.017 

0.043 

0.076 

0.123 

±0.035 

±0.063 

±0.096 

±0.143 

1 

±0.029 

±0.043 

±0.062 

±0.097 

O-O32 

0.054 

0.081 

0.125 

-0.038 

-0.073 

-0.124 

-0.219 

+0-043 

+0.082 

+0.135 

+0.222 

±0.024 

±0.058 

±0.112 

±0.217 

0.013 

0.031 

0.057 

0.100 

-0.006 

-0.024 

-0.055 

-0.117 

±0.014 

±0.029 

±0.050 

±0.091 

+0.006 

+0.008 

+0.010 

+0.017 

-0,018 

-0.044 

-0.086 

-0.167 

0.009 

0.024 

0.053 

0.122 

E  2.4 

S  FOR  \  =  0.625 


r  ■  ..../.o 

X  /c  =  0.0  X  /c  =  ±0.125 

P  P 


\  fx 
n 

0.90 

0.94 

0.97 

- \ 

1.00 

0.90 

0.94 

0.97 

- \ 

1.00 

0 

1.003 

1.072 

1.123 

1.17.5 

1.007 

1.077 

1.129 

1.183 

1 

0 

0 

0 

0 

±0.340 

±0.369 

±0.392 

±0.4l4 

2 

0.659 

0.747 

0.820 

C.899 

0.567 

0.645 

0.710 

0.781 

3 

0 

0 

C 

0 

+0.288 

+0.343 

+0.391 

+0.445 

4 

-0.284 

-0.354 

-0.417 

-0.491 

-0.149 

-0.186 

-0.219 

-0.260 

5 

0 

0 

0 

0 

±0.188 

±0.248 

±0.306 

±0.378 

6 

0.148 

0.2C4 

0.260 

0.331 

0.010 

0.012 

0.014 

0.018 

7 

0 

0 

0 

0 

+0.106 

+0.156 

+0.209 

+0.279 

8 

-0.086 

-0.131 

-0.181 

-0.249 

0.035 

0.055 

0.078 

0.109 

9 

0 

0 

0 

0 

±0.050 

±0.081 

±0.117 

±0.169 

10 

0.053 

0.090 

0.134 

0.200 

-0.041 

-0.071 

-0.108 

-0.163 

11 

0 

0 

0 

0 

+0.015 

+0.027 

+0.041 

+0.064 

12 

-0.034 

-0.064 

-0.103 

-0.166 

0.032 

0.061 

0.101 

0.165 

TABLE 
Si 


COEFFICIENTS 


±0.673 

±0.737 

±0.786 

±0.835 

0.297 

0.341 

0.380 

0.424 

+0.401 

+0.485 

+0.560 

+0.647 

0.129 

0.170 

0.207 

0.250 

±0 . 096 

±0 . 126 

±0.157 

±0.197 

-0.132 

-0.190 

-0.251 

-0.332 

±0.045 

±0.073 

±0.102 

±0.143 

0.040 

0.061 

0.087 

0.125 

+0.056 

+0.097 

+0.146 

+0.222 

0.019 

0.038 

0.062 

0.100 

±0.019 

±0.034 

±0.055 

±0,091 

-0.026 

-0.055 

-0.096 

-0 .166 

:.E  2.5 

TS  FOR  X  = 

0.75 

1.011 

1.101 

1.171 

1 

.242 

±0.967 

±1.083 

±1.175 

±1 

.270 

-0.109 

-0.138 

-0.158 

-0, 

.175 

+0.209 

+0.253 

+0.298 

+0 

.358 

0.222 

0.299 

0.377 

0 

.480 

+0.126 

+0.191 

+0.261 

+0 

.356 

0.026 

0.053 

0.083 

0 

.123 

±0.034 

±0.048 

±C.066 

±0 

.097 

-0.049 

-0.087 

-0.136 

-0 

.219 

±0.035 

±0.072 

±0.125 

±0 

.217 

-0.011 

-0.031 

-0.062 

-0 

.117 

+0.007 

+0.009 

+0.011 

+0 

.017 
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